Photoionization effects in streamer discharges 

Manuel Arrayas 1 , Marco A. Fontelos 2 and Jose L. Trueba 1 
1 Departamento de Fisica Aplicada, Universidad Rey Juan Carlos, 
Tulipdn s/n, 28933 Mostoles, Madrid, Spain and 
2 Instituto de Matemdticas y Fisica Fundamental, 
Consejo Superior de Investigaciones Cientificas, 
C/ Serrano 123, 28006 Madrid, Spain 
(Dated: February 2, 2008) 

Abstract 

In this paper we study the effects of photoionization processes on the propagation of both 
negative and positive streamer discharges. We show that negative fronts accelerate in the presence 
of photoionization events. The appearance and propagation of positive streamers travelling with 
constant velocity is explained as the result of the combined effects of photoionization and electron 
diffusion. The photoionization range plays an important role for the selection of the velocity of the 
streamer as we show in this work. 

PACS numbers: 52.80.-s, 94.05.-a, 51.50.+V 
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I. INTRODUCTION 

Since Raether |j] used cloud chamber photographs to study the creation and propaga- 
tion of streamer discharges there has been considerable effort to understand the underlying 
processes driving them. A streamer discharge is considered to be a plasma channel which 
propagates in a gas. The discharge propagates by ionizing the medium in front of its charged 
head due to a strong field induced by the head itself. This kind of discharges produces sharp 
ionization waves that propagate into a non-ionized gas, leaving a non-equilibrium plasma 
behind. 

Raether himself realized that Townsend's mechanism which takes into account the cre- 
ation of extra charge by impact ionization Q was not enough to explain the velocity of 
propagation of a streamer discharge. He pointed to photoionization as the process which 
enhances the propagation of the streamer. Due to the recombination of positive ions and 
electrons, the head of the discharge is a strong source of high energetic photons. Photons, 
emitted by the atoms that previous collisions have excited, initiate secondary avalanches in 
the vicinity of the head which move driven by the local electric field increasing the velocity 
of propagation of the front. 

In this paper we study the role played by photoionization in the propagation of both 
negative and positive streamers. We take a model widely used in numerical simulations 
and find an effective simplified model. We discuss how this simplified model retains all the 
physics of streamer discharges including photoionization. The photoionization is modelled 
as a nonlocal source term. We take the case of air and consider optical emissions from N 2 
and molecules. Then we consider the sole role of photoionization in negative planar 
shock fronts. Finally we analyse the case of positive planar fronts and propose a mechanism 
for their formation and propagation. We end with an analysis of results and conclusions. 

II. MODEL FOR A STREAMER DISCHARGE 

Here we consider a fluid description of a low-ionized plasma based on kinetic theory. 
The balance equation for the particle density of electrons N e is the lowest moment of the 
Boltzmann equation, 




(1) 
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where R is the position vector, r is time, Vr is the gradient in configuration space, U e (R, r) 
is the average (fluid) velocity of electrons and S e is the source term, i.e. the net creation rate 
of electrons per unit volume as a result of collisions. It is convenient to define the electron 
current density J e (R, r) as 

J e (R,r) = iV e (R,r)U e (R,r), (2) 

so that the balance equation can also be written as 

BN 

^ + V R -J e = S e . (3) 

The same procedure can be done, in principle, for positive (N p ) and negative (N n ) ion 
densities to give 

ON 

^ + V R -J P = S P , (4) 
^ + V R -J n = S n , (5) 

where J Pin are the current densities of positive and negative ions, respectively, and S PtTl 
are source terms. Conservation of charge has to be imposed in all processes, so that the 
condition 

S p — S e + S n , (6) 

holds for the source terms. Some physical approximations can now be done in order to 
simplify the balance equations (p3|)-(p)f . The first one is to assume that the electron current 
J e is approximated as the sum of a drift (electric force) and a diffusion term 

J e = -fx e SN e - D e V R iV e , (7) 

where £ is the total electric field (the sum of the external electric field applied to initiate the 
propagation of a ionization wave and the electric field created by the local point charges) and 
H e and D e are the mobility and diffusion coefficient of the electrons. Note that, as the initial 
charge density is low and there is no applied magnetic field, the magnetic effects in equation 
(JZJ) are neglected. Concerning the diffusion coefficient, in the case of equilibrium, the kinetic 
theory of gases links diffusion to mobility through Einstein's relation D e / /j, e = kT/e. With 
respect to positive and negative ions, on time-scales of interest for the case of streamer 
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discharges, the ion currents can be neglected because they are more than two orders of 
magnitude smaller than the electron ones Q], so we will take 



Jp = Jn = 0. (8) 

Consider now the processes that give rise to the source terms S e , p ,n- 

1. The first of these processes is the creation of free electrons by impact ionization: an 
electron is accelerated in a strong local field, collides with a neutral molecule and 
ionizes it. The result is the generation of new free electrons and a positive ion. The 
ionization rate is given by 

SI = S l p = UiN e , (9) 

where the ion production rate z/, depends on the local electric field, the density of the 
neutral particles of the gas and their effective ionization cross sections. 

2. The second possible process is attachment: when an electron collides with a neutral 
gas atom or molecule, it may become attached, forming a negative ion. This process 
depends on the energy of the electron and the nature of the gas [9]. The attachment 
rate can be written as 

S a n = -S a e = u a N ej (10) 

where v a is the attachment rate coefficient. Note that the creation of negative ions 
due to these processes reduces the number of free electrons, so S® is negative. 

3. There are also two possible kinds of recombination processes: a free electron with a 
positive ion and a negative ion with a positive ion. The recombination rate is 

Sf = S e p p = -u ep N e N p , (11) 

for electron-positive ion recombination, and 

S;p = S r : p = -u np N n N p , (12) 

for positive ion-negative ion recombination, v ep and v np being the recombination coef- 
ficients respectively. 
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4. Finally, we can include photoionization: photons created by recombination or scatter- 
ing processes can interact with a neutral atom or molecule, producing a free electron 
and a positive ion. Models for the creation rate of electron-positive ion pairs due to 
photoionization are non-local. This rate will be here denoted by 

Sf = S ph = S ph . (13) 

Taking into account the expressions (J7J) and (jS)) for the current densities, and the equations 
(J§J)- (jl3|) for the source terms, we obtain a deterministic model for the evolution of the 
streamer discharge, 

dN P . 



dr 



= V R • (fi e SN e + D e V R N e ) + UiN e 

- v a N e -v ep N e N p + S ph , (14) 

O TV T 

_JL = ViNe _ Uep N e N p - v np N n N p + S ph , (15) 

_ZL = UaNe _ y np N n N p . (16) 

In order for the model to be complete, it is necessary to give expressions for the source 
coefficients u, the electron mobility /i e , the diffusion coefficient D e and the photoionization 
source term S ph . Finally, we have to impose equations for the evolution of the electric field 
S. This evolution of the electric field is given by Poisson's equation, 

V R ■ S = - (N p - N n - N e ) , (17) 

where e is the absolute value of the electron charge, e is the permittivity of the gas, and we 
are assuming that the absolute value of the charge of positive and negative ions is e. Note 
that the coupling between the space charges and the electric field in the model makes the 
problem nonlinear. The model given by ()14j) . (|15p. and (|16jh together with ([17)1 has been 
studied numerically in the literature [l(| . There are other works where the electrical current 
due to ions (jH)) is taken into account although not photoionization jisl ]. 



III. A SIMPLIFIED MODEL 

In this section we will simplify the model given by equations ()14|) - ([T6|) . In order to be 
specific and fix ideas we shall consider the case of air. In [n]], some data are presented for 
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the ionization coefficients and the photoionization source term. Using these data we shall 
see that one can neglect the quadratic terms involving the coefficients v ep and v np since they 
are about two orders of magnitude smaller than z^. The same can be said about the terms 
involving the coefficient v a . First we write equations (|14j) - (jl6j) as 

dN 

= V R - (ve£N e + D e V R N e ) 

+ ( Ui -u a -u ep N p )N e + S ph , (18) 



dr - = & - v ep N p ) N e - v np N n N p + S ph , (19) 
dN 

= v a N e - v np N n N p . (20) 

In these equations, and using the data in (Figure 1 and Table 2), the term Ui is of 
the order of 10 10 s _1 for large electric fields, v a is about 10 8 s _1 , and v ep and v np are about 
10~ 13 m 3 ■ s _1 . Moreover, N p is of the same order of N e . Then, in equation (12 OH . in the 
stationary regime when the particle densities reach the saturation values, one has N n ~ 
v al v np ~ 10 21 m~ 3 . So that, it follows from equation ()19|) that, in the stationary regime, the 
term u np N n N p ~ 10 8 A^ p is two orders of magnitude smaller than the term z^iV e ~ 10 10 iv" e . 
Hence the terms v a N e and v np N n N p can safely be neglected. The model then reads 

dN 

— ± = V R ■ (fi e SN e + D e V n N e ) 
or 

+ (u l -u ep N p )N e + S ph , (21) 
^ = (u i -u ep N p )N e + S ph . (22) 

In order to neglect the term v ep N e N p by comparison with the term UiN e , it is necessary than 
N p (and then iV e ) satisfies N p Vijv ep ~ 10 23 m -3 . To see that it is the case, we use the 
Poisson equation (|17j) to write equation (J21|) . without the term v ep N e N p , as 

' ' - fi e £-V R N e -D e V 2 n N e 



Or 



(vi + li e ^-(N p -N e )^N e + S ph . (23) 



From this expression, looking at its RHS, we can see that, while S ph has small effect and 
the total populations of both ions and electrons, N e can grow only up to a saturation value 
at which z/j + /i e — (N p — N e ) = 0, i.e. 



■so 

v 



N e -N p < — 10 20 m 3 , (24) 

fJL e e/£ 



at all times. Therefore neither N p nor N e reach values close to 10 23 m 3 , and all the assump- 
tions which led to neglect i> ep N e N p are justified. Our simplified model will be 

C\ TV T 

-^=V R - (fi e £N e + D e VnN e ) + v,N e + S ph , (25) 
^ = Vi N e + S* h . (26) 

Let us remark that the orders of magnitude deduced for N e and N p coincide with those 
found in full numerical simulations by Liu and Pasko 1(3 ] . 

IV. THE PHOTOIONIZATION TERM 



In this section we will write down an explicit form of the photoionization source term. 
In our study on the effects of photoionization on the evolution of streamers in air we con- 
sider that only optical emissions from N2 and molecules can ionize O2 molecules. The 
photoionization rate, due to the fact that the number of photons emitted is physically pro- 
portional to the number of ions produced by impact ionization, is written as the following 
nonlocal source term jloL lri|. 

S ph (R) = S J Vi(R')N e (B!) K ph (\R - R'|) d 3 R', (27) 

where So is given by 

s ° = r^(-) , < -1 v < 28 > 

max X'min I 

In this expression, p q is the quenching pressure of the single states of N 2 , p is the gas pressure, 
£ is the average photoionization efficiency in the interval of radiation frequencies relevant 
to the problem, v* is the effective excitation coefficient for N2 state transitions from which 
the ionization radiation comes out (we take v*/vi to be a constant), and Xmin and Xmax are, 
respectively, the minimum and maximum absorption cross sections _of O2 in the relevant 
radiation frequency interval. The kernel i^ p /j(|R — R'|) is written as 



12 



Krk(R) = exp (-X.fi) - 3 exp(-, 2fl) i (2g) 

in which Xi — XminPo-z and X2 = XmaxPon so that xi < X2- F° r the ionization coefficient z/j, 
we take the phenomenological approximation given by Townsend j^, 

Ui = ii e \£\a exp( -—f) , (30) 



where fi e is the electron mobility, «o is the inverse of ionization length, and So is the char- 
acteristic impact ionization electric field. Note also that /J, e \£\ is the drift velocity of elec- 
trons. Townsend approximation provides some physical scales and intrinsic parameters of 
the model. It is then convenient to reduce the equations to dimensionless form. Natural 
units are given by the ionization length R = a^ 1 , the characteristic impact ionization field 
Sq, and the electron mobility /i e , which lead to the velocity scale Uq = [i e So, and the time 
scale T = Rq/Uq. We introduce the dimensionless variables r = H/Rq, t = t/tq, the 
dimensionless field E = B/Sq, the dimensionless electron and positive ion particle densi- 
ties n e = N e /N and n p = N p /N with N = EoSo/(eRo), and the dimensionless diffusion 
constant D = D e /(R U ). The dimensionless model reads then, 

= V • (n e E + DVn e ) + n e |E|e~ 1/|E| + S, (31) 

^ = n e |E|e- 1 /l E l+5. (32) 

where S is the dimensionless photoionization source term, 

S(r) = S J ^(rOlE^OIe-W^lirdr-r'DdV, (33) 

g = j— gg-£(^ V ? \ - v (34) 

max I X'inin J 

K ( r j = exp (-(xi/a )r) - exp (-(xi/aojr) ^ 

In this paper, we restrict ourselves to a planar geometry, in which the evolution of the 
ionization front is along the z-axis. In this case, the photoionization source term can be 
written as 

S(z)=S J dz' n e (z' \t)\E(z' \t)\e~ mz '^ I(\z - z'\), (36) 

where 

/oo roc | 
d y' / dx> i ITi /2 i 7 AM 
y_oo [x 1 + y' 2 + {z- z'yy/ 2 

x ^ e (-(xi/ao)V*' 2 +w' 2 +( z -*') 2 ) _ e (-(X2M)) V /;c ' 2 +f' 2 + (^ z ') 2 )j _ (37) 

Changing to cylindrical coordinates, and integrating in the polar angle, equation (|37|) results 

exp (-(xi/«o)\A J + (z - z') 2 ) - exp (-(x2/«o)\A" 2 + (z - z') 2 ) 

(r 2 + {z- z') 2 ) 3 / 2 

(38) 



and 



Also, 



in 

POO 

I(\z-z'\) = 2tt / rdr 
Jo 



We can define s — \z — z'\ and w = -y/r 2 + s 2 . Then, 



I(s) = 2n r dw «P(-(Xi/°o)">)-«p(-(X2/ao)«>) (3g) 

Js W 

Denning the quantities 

= 2nS = \J£- i ( V A -j—L (40) 

max Xmin ) 

and 

k(s) = (41) 
we can write the dimensionless photoionization term in the planar case as 

S(z) = Vo J dz'n e (z',t)\E(z , ,t)\e- imz '' t)l k(z-z'), (42) 

where 

k(s) = r dx ^(-Xix)-^V(-X2x)_ (43) 

Js/ ao a X 2 

The function k(s) cannot be computed explicitly in terms of elementary functions, but its 
asymptotic behaviour can be calculated. For s — >• oo, we have 

e -0»Ci/ao)s e -(X'2/ao)s 

k ^ - ( I \ 2 - ( I \ 2 ' ( 44 ) 

and for s — > 0, it is 

k(s) ~ Xl ~ X2 In s + const. (45) 

In the numerical computations, we will approximate the function k(s) by functions with 
the same behaviour at infinity and zero as the ones shown in equations (jUj) and ([45)1 . The 
simulations show that the result is insensitive to the details of these approximations and 
they only depend on the behaviour at zero and infinity. In fact, we will use a kernel such 
that it is equal to (|43|) for s < 1 and it is equal to (}4*4"|) for s > 1. The constant in equation 
()4*5|l will be chosen in such way that k(s) is continuous at s = 1. 

Following and [l8|, we will take for the simulations = 0.1, p g = 30Torr, 

Xi = 0.035 Torr _1 cm _1 p 02 , X2 = 2 Torr _1 cm~ 1 p 02 . We will assume the partial pressure 
of the oxygen in air is given by po 2 — IP, where p is the total pressure and 7 a pure 
number between zero and one. For the inverse ionization length ao, we will take the value 



for nitrogen, that depends on pressure 
coefficient jl3j], we take D e = 0.1m 2 /s. 



as «o = 5.8Torr 1 cm 1 p. For the diffusion 



Using these values it turns out, 



«. = 0.37^, (46) 



with p expressed in Torr, and 



cxp(-0.006 7£) _ cxp(-0.34 7s) 
k ( s ) = J (0.006 7 )s 2 (0.34 7 ) S 2 ' S > 1, 

O^lns I cx P(-°- 006 ^ cxp(-0.34 7 ) <-, V 7 

U. ^71116 -f (0.006 7) (0.34 7 ) ' 6 - X > 

V. PHOTOIONIZATION WITHOUT DIFFUSION: ACCELERATION OF NEGA- 
TIVE FRONTS 

We consider the case in which a divergence-free electric field Eo = —Equ z is set along the 
z-axis, so that electrons move towards the positive z-axis. Then we take the electric field 
as E = —Eu z , E being its modulus, so that, in the case in which the diffusion coefficient is 
D — 0, the model can be written as 

^ = -^{n e E)+n e Ee^ E + S } (48) 

^ = n e Ee- x ' E + S, (49) 

dE . . 

n p -n e = - — . (50) 



Now, following the approach presented in 



16, 



17(, we introduce the shielding factor u(z,t) 



as 



in terms of which, 



u{z,t) =e-ti neiz ' t,)dt ', (51) 



1 du 

n e = - — , 52 
u at 

1 du dE u 

n * = ~um ~ ~dz~' (53) 

E = E u, (54) 



and hence 



S(z) = ipo I dz' n e {z')E^z')u{z')e- xjE ^ z ' )u ^ k(z - z') = -ip ^_ J dz' G(u(z')) k(z - z'), 

(55) 
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where 

G{u) = -f d Ul E e~ 1/EoUl . (56) 
u I — 1 1 — 
In order to deduce an equation for the shielding factor u, we follow the steps of jig, |l7[ and 

obtain a Burgers equation with non-local source 

Ou du f 

— + E u— = -uripo + uG(u) + (f u / G{u{z')) k(z - z'), (57) 
u(z,0) = 1, (58) 

where n p o is the initial positive ion density. Our method of solution of the above system is 
by integration along characteristics; i. e. we solve the following system of ODE's 
dz 

* = E " % (59) 

— = — n p0 u + uG(u) + cpou / dz G(u(z'))k(z — z'). (60) 

We use this formulation in terms of characteristics in order to give a numerical algorithm 
and study the effect of photoionization on the propagation of negative planar fronts. We 
discretize the spatial variable z into N segments separated by the points Zo, z%, ... zn and 
follow the evolution in time of each of them by solving (j5T?j) and (I60|) . The integral term in 
(160 j) is discretized in the following form 

„ 2V-1 

/ dz'G(u(z , ))k(z-z , )~J2 G (< z ^ k ( z ^- z ^) (^i+iW -%(*))■ ( 61 ) 

J j=0 

In our first numerical experiment, we choose as initial data a Gaussian distribution of 
charge. We take E = 1 and the pressure p = 750 Torr. In FigCJwe can see the evolution of 
the initial negative charge distribution when the photoionization term is neglected. It can 
be seen that electrons move in the direction of increasing z where the anode is situated. A 
negate front is deve.oped at the rig ht of the initial distribution Q. The e.ecttons at the 
left side of the initial distribution move also following the electric field, until they reach the 
main body of the plasma where the electric field is screened. Then they stop there (around 
z = 2 in FigU). When the photoionization term is included, the profiles change. In Fig|21 
the same numerical experiment is carried out, with the inverse of photoionization range 
7 = 0.21, which corresponds to the normal conditions of air in the atmosphere. 

We can track the motion of the negative front by looking at the time evolution of the 
point z*(t) at which the electron density has a given value. In Fig|21 we compare the graphs 
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FIG. 1: Electron density n e profiles without photoionization. The electrons move to the right 
following the polarity of the electric field. A negative planar front is developed. 

4.5i 1 1 1 1 1 1 



FIG. 2: Electron density n e profiles with photoionization, at normal pressure p = 750Torr and 
7 = 0.21. A negative front is moving towards the anode at the right and electrons start getting 
accumulated at the zero electric field plasma zone. 

of z*(t) with and without photoionization for a level of n e = 0.1. As we can see, the effect of 
photoionization is an acceleration of the negative front which reaches a higher though still 
constant velocity. This fact holds, after our observations, when one considers kernels k(s) 
which decay exponentially fast at infinity. 

Finally, it is interesting to observe the behaviour of the density n e in the direction opposed 




z 
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3 1 — i , , 1 1 , 1 , , — 

0.4 0.6 0.8 1 1.2 1.4 1.6 1.8 2 
t 

FIG. 3: The evolution of a point z* of the negative front at which the electron density has the 
value n e = 0.1. When photoionization range 1/7 is increased, the front moves faster. The line 
without label belongs to the case where photoionization is neglected. 

to the propagation of the negative front (the left part of the initial distribution). This will 
be called from now on "the positive front" . We can observe in FigEl an effect consisting in 
the accumulation of electrons in a small region of space in the positive front. This fact is 
easy to understand by considering the production of electrons away from the positive front 
which are drifted towards the positive front following the electric field. In the positive front, 
electrons and positive ions are balanced and hence the net electric field cancels. Therefore 
electrons cannot proceed any further beyond the positive front and they accumulate there. 
This is an effect purely associated to photoionization which cannot be explained by invoking 
any different effect. Unless there is some mechanism allowing the electrons to spread out 
once they accumulate at the positive front, their density will grow indefinitely and eventually 
will blow up. We will see in next section that this mechanism is diffusion and the net effect 
of photoionization and diffusion is the appearance of travelling waves moving towards the 
cathode, i.e. positive ionization fronts. 
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FIG. 4: Electron density n e profiles, at normal pressure p = 750Torr, photoionization parameter 
7 = 0.25 and diffusion D = 0.57 in dimensionless units. A negative front is moving towards the 
anode at the right and a positive fronts towards the cathode at the left. 

VI. PHOTOIONIZATION WITH DIFFUSION: POSITIVE IONIZATION 
FRONTS 

In this section we study in one space dimension the combined effect of photoionization 
and diffusion on the propagation of positive fronts. The system of equations we study is 
therefore 

^--l(n e E- D ^) + n e Ee^ + S, (62) 
^ = n e Ee'" E + S, (63) 

n p~ n e = — qJ, (64) 

where S is the photoionization source term and is written as in equation (|55|). 

The main difference in our approach to this problem with respect to the problem without 
diffusion is that now an integration along characteristics does not lead to simplifications due 
to the presence of the second derivatives associated with diffusion. Instead we will use the 
method of finite differences. 

In FigEJ we represent the profiles for n e with D = 0.57, p = 750Torr and 7 = 0.25. We 
have used an initial charge distribution which has a maximum at z = 10. When it evolves, it 
can be observed a negative planar front developing. The propagation of the negative front is 
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FIG. 5: The evolution of points z* of the negative and positive fronts at which the electron density 
has the value n e = 0.02. The increasing values are for the negative front and the decreasing ones 
are for the positive. When photoionization range I/7 is increased, the fronts move faster. Triangles 
A are for 7 = 0.9, diamonds for 7 = 0.25 and squares □ for 7 = 0.1. 

almost identical with or without diffusion when photoionization is present. However there is 
now a positive front moving towards the cathode. The positive front moves with a constant 
velocity which is smaller than the velocity of the negative front. In Fig IS] we have plotted the 
position z* of a point of the negative front and of the positive front which has the particular 
value of the electron density n e = 0.02. The parameters are the same as in Fig0J but for 
three different values of 7. For the parameter values chosen above, we have computed the 
ratio between the velocities of positive and negative fronts: c pos /c neg = 0.34 for 7 = 0.9, 
Cpos/cneg = 0.68 for 7 = 0.25 and c pos /c neg = 0.86 for 7 = 0.1. The ratio grows when the 
photoionization range I/7 increases and the velocities for negative and positive fronts tend 
to increase and get closer to each other. 

The propagation of positive fronts as travelling waves results from the combined action 
of photoionization and diffusion. This is in contrast with the propagation mechanism for 
negative fronts, which are also travelling waves but they result from a combination of impact 
ionization and convection by the electric field. In the latter case, diffusion and photoion- 
ization only affect the negative fronts by changing their velocity and their shape. All this 
conclusions are rather insensitive to the detailed form of the kernel k(s) (see formula (|41j) ) 
provided it decays exponentially fast at infinity, and hence our conclusions hold with a high 
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degree of generality. 



VII. CONCLUSIONS 

In this paper we have studied the effect of photoionization in streamer discharges. We 
have deduced a minimal model including photoionization and studied with this model the 
propagation of both positive and negative fronts in the planar case. We have found the ap- 
pearance of travelling waves which accelerate when the photoionization range increases. For 
negative fronts we have studied the effect of photoionization both when electronic diffusion 
is neglected and included. For positive fronts, electronic diffusion has to be taken into ac- 
count and we have shown how photoionization plays the crucial role pointed by Raether on 
increasing the velocity of propagation. The control parameter is the photoionization range, 
i.e. the typical distance at which photons are able to ionize the media. Physically in air, 
this parameter depends on the amount of oxygen and nitrogen present. It is interesting to 
point out that for real discharges in the atmosphere, this parameter varies with the altitude. 



[1] H. Raether, Die Entwicklung der Elektronenlawine in den Funkenkanal, Z. Phys. 112, 464-489 
(1939). 

[2] L. B. Loeb, The problem of the mechanism of static spark discharge, Rev. Mod. Phys. 8, 
267-293 (1936). 

[3] M. Arrayas and J. L. Trueba, Investigations of Pre-Breakdown Phenomena: Streamer Dis- 
charges, Cont. Phys. 46, 265-276 (2005). 
[4] Y. P. Raizer, Gas Discharge Physics (Springer, Berlin 1991). 

[5] L. B. Loeb and J. M. Meek, The mechanism of the electric spark, Clarendon Press, Oxford, 
1941. 

[6] V. P. Pasko, M. A. Stanley, J. D. Mathews, U. S. Inan, and T. G. Wood, Electrical discharge 
from a thundercloud top to the lower ionosphere, Nature 416, 152-154 (2002). 

[7] M. Arrayas, U. Ebert, and W. Hundsdorfer, Spontaneous branching of anode-directed stream- 
ers between planar electrodes, Phys. Rev. Lett. 88, 174502 (2002). 

[8] M. Arrayas, M. A. Fontelos, J. L. Trueba, Mechanism of branching in negative ionization 



16 



fronts, Phys. Rev. Lett. 95, 165001 (2005). 
[9] S. K. Dhali and A. P. Pal, Numerical simulation of streamers in SFq, J. Appl. Phys. 63, 
1355-1362 (1988). 

[10] N. Liu and V. P. Pasko, Effects of photoionization on propagation and branching of positive 
and negative streamers in sprites, J. Geophys. Res. 109, A04301 (2004). 

[11] G. V. Naidis, On photoionization produced by discharges on air, Plasma Surces Set. Technol. 
15, 253-255 (2006). 

[12] M. B. Zhelezniak, A. Kh. Mnatsakanian, S. V. Sizykh, Photoionization of nitrogen and oxygen 
mixtures by radiation from a gas discharge, High Temperature, 20, 357-362 (1982). 

[13] P. A. Vitello, B. M. Penetrante, and J. N. Bardsley, Simulation of negative-streamer dynamics 
in nitrogen, Phys. Rev. E 49, 5574-5598 (1994). 

[14] U. Ebert, W. van Saarloos, and C. Caroli, Streamer propagation as a pattern formation 
problem: Planar fronts, Phys. Rev. Lett. 77, 4178-4181 (1996), and ibid., Propagation and 
structure of planar streamer fronts, Phys. Rev. E 55, 1530-1549 (1997). 

[15] M. Arrayas, On negative streamers: A deterministic approach, Am. J. Phys. 72(10), 1283- 
1289 (2004). 

[16] M. Arrayas, M. A. Fontelos, J. L. Trueba, Ionization fronts in negative corona discharges, 

Phys. Rev. E 71, 037401 (2005). 
[17] M. Arrayas, M. A. Fontelos, J. L. Trueba, Power laws and self-similar behaviour in negative 

ionization fronts, J. Phys. A: Math. Gen. 39, 1-18 (2006). 
[18] A. A. Kulikovsky, The role of photoionization in positive streamer dynamics, J. Phys. D: 

Appl. Phys. 30, 1514-1524 (2000). 



17 



